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Abstract. The expected chirp mass distribution of ob- 
served events for future gravitational wave detectors is 
extensively investigated in the presence of an exotic fluid 
component with an arbitrary equation of state, — 1 < ui x = 
Px/px < 0, i.e., the so-called dark energy component. The 
results for a flat model dominated by a dark energy are 
compared to those for the standard flat model dominated 
by cold dark matter. It is found that for a flat universe 
the chirp mass distribution shows a sensitive dependence 
on lo x , which may provide an independent and robust con- 
straint on the cosmic equation of state. 
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1. Introduction 

Upcoming interferometric detectors of gravitational waves 
such as LIGO (Abramovici et al. 1992), VIRGO (Caron 
et al. 1997), GEO (Danzmann et al. 1997) and TAMA 
(Kuroda 1997) will produce a flood of data when they 
come online in a few years. Advanced detectors in 
LIGO/VIRGO and LCGT (Kuroda et al. 1999) will even 
be able to see neutron star binaries (NS/NS) out to cos- 
mological redshifts. Because the accuracy in the measure- 
ment of the signal strength can be better than 10%, and 
the accuracy in the measurement of the chirp mass (which 
characterizes the frequency sweep of a binary inspiral) 
can be better than 0.1%(Finn & Chernoff 1993; Cutler 
& Flanagan 1994), observations of NS/NS inspiral events 
can be used to determine the universe's Hubble constant, 
deceleration parameter and cosmological constant (Schutz 
1986, 1989; Markovic 1993; Chernoff & Finn 1993; Finn 
1996; Wang & Turner 1997). In particular, Finn (1996) 
and Wang & Turner (1997) explicitly pointed out that 
the chirp mass spectrum of observed neutron star binary 
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inspiral events can be quite sensitive to cosmological con- 
stant. 

In this paper, we investigate the possibility of using the 
future gravitational wave data to constrain dark energy, an 
exotic fluid component in the universe with an arbitrary 
equation of state, —1 < u> x = p x /p x < 0, whose existence 
has been well motivated by various astronomical observa- 
tions and theoretical arguments (refer to a recent review, 
Turner 1999). While the cosmological constant (vacuum 
energy) can apparently succeed in explaining the acceler- 
ating expansion of the universe which has stirred cosmol- 
ogists after observations and measurements of distant la 
type supernovae (SNela) (Perlmutter et al. 1999; Riess et 
al. 1998), it suffers from the difficulty in understanding of 
the observed value in the framework of modern quantum 
field theory (Weinberg 1989; Carroll et al. 1992) and the 
"coincidence problem" , the issue of explaining the initial 
conditions necessary to yield the near-coincidence of the 
densities of matter and the cosmological constant compo- 
nent today. In this case, a dark energy component with 
generally negative pressure has been invoked. Examples 
of dark energy include the evolving scalar field (Caldwell 
et al. 1998; Zlatev et al. 1999; Steinhardt et al. 1999), the 
smooth noninteracting component (Turner & White 1997; 
Chiba et al. 1997) and the frustrated network of topologi- 
cal defects in which uj x — — §, with n being the dimension 
of the defect (Spergel & Pen 1997). Some observational 
constraints for these models have been extensively ana- 
lyzed in the literature, including the redshift-luminosity 
distance relation based on distant SNela (Garnavich et al. 
1998), gravitational lensing statistics (Cooray & Huterer 
1999; Zhu & Cao 1999; Zhu 2000), Einstein ring system 
(Futamase & Hamana 1999), the anisotropy of the Cos- 
mic Microwave Background Radiation (CMBR) (Huey et 
al. 1999), the angular size-redshift relation of distant ex- 
tragalatic sources (Lima & Alcaniz 2000a), the age of 
high-rcdshift galaxies (Lima & Alcaniz 2000b), and the 
combinations of SNela measurements with either gravi- 
tational lensing (Waga & Miceli 1999) or CMBR (Efs- 
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tathiou 1999) or large-scale structure (Perlmutter et al. 
1999) or the time variations of the redshift of quasars 
(Nakamura & Chiba 1999). A combined maximum like- 
lihood analysis has suggested a range of equation-of-state 
for quintessence, —1 < w x < —0.6 (see Wang et al. 2000 
for an elegant summary). 

Because the determination of the amount and nature of 
dark energy is emerging as one of the most important chal- 
lenges in cosmology, it is worthwhile to investigate other 
observational aspects of its equation of state. In this work, 
from calculating the expected chirp mass spectrum of ob- 
served NS/NS inspiral events for future gravitational wave 
detectors in the presence of a dark energy component, we 
find a sensitive dependence of chirp mass distribution on 
ui x for a flat universe, which could provide an independent 
and robust constraint on the cosmic equation of state. 
For calculations of this work, our notation and approach 
mainly follow those described by Finn (1996). 

2. The chirp mass spectrum of NS/NS inspiral 
events 

The NS/NS inspirals are the targets for future gravita- 
tional wave detectors, whose merger rate at redshift z per 
unit observer time interval per unit physical volume is, 
n = hoF(z) (1 + z) 2 , where F(z) describes the NS/NS 
merge rate evolution, and no(— 10 _7 /iMpc _3 yr~ 1 ) is the 
local NS/NS merger rate per unit volume (Phinney 1991). 
In the case of no proper evolution F(z) = 1. The evolu- 
tion of the NS/NS merger rate has been studied by several 
authors (e.g., Tutukov & Yungclson 1993; Lipunov et al. 
1995). Because recent observations of field galaxies out to 
redshift z ~ 5 (see e.g. Madau et al. 1998; Stcidcl et al. 
1999) have provided a means of modeling the evolution of 
the star formation rate (SFR), one can now compute the 
birth and merger rates of binary systems by taking into 
account the time-delay. Using population synthesis model, 
Schneider et al. (2000) have calculated birth and merger 
rate evolution of various kinds of binaries. Our best fit 
to their results for NS/NS merger rate from hierarchical 
scenario (Fig. 8 of Schneider et al. 2000) takes 

F(z) = l + 11.5z 2 -13.6z 3 +5.37z 4 -0.70z 5 ze(0,3).(l) C e (x) = I 



the detector bandwidth, while tq gives an overall sense of 
the depth to which the detector can reach. The signal- 
to-noise ratio of an inspiraling binary with a single inter- 
feromctric detector depends on the relative orientation of 
the source, which is confined to the angular orientation 
function 6. From observations of binary inspiral in a sin- 
gle interferometer we can measure p and M. but not 0. 
Fortunately, even though G cannot be measured, its prob- 
ability distribution can be approximated by (Finn 1996) 
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which is related to P p (p\A4o,Ai), the probability distri- 
bution of signal-to-noise ratio of events with chirp mass 
M, through 
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where Z = M.j '.Mo — 1- Then, in a homogeneous and 
isotropic universe we get the rate at which binary inspiral 
signals corresponding to chirp mass M. and p greater than 
Po are detected 
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where Da is the angular diameter distance from the de- 
tector to the source, Cq(x) is the probability that a given 
detector receives a binary inspiral with chirp mass M. and 
signal-to- noise ratio greater than p§: 



(1 + x) (4 - a;) 4 /256, 
0, 



if < x < 4, 
otherwise. 
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For a NS/NS inspiral source located at redshift z, 
the detectors measure its chirp mass, M. = .Mo (1 + z ) 
where the intrinsic chirp mass, .Mo = /i 3 ^ 5 M 2 / 5 , with 
fi and M being the binary's reduced and total mass 
respectively. Therefore, the chirp mass spectrum repre- 
sents directly the source redshift distribution if .Mo is 
constant. Because we are considering NS/NS binaries, 
Mo = 1.19 M Q will be used below, corresponding to a 
typical neutron star mass m = 1.37 M Q . The signal-to- 
noise ratio p of an event observed by a given detector is 
p = 86 (r /D L ) (M/1.2M f 6 C(/max), where D L is the 
luminosity distance from the detector to the source. £(~ 1 
for our case) describes the overlap of the signal power with 



An accumulation of observed NS/NS binary inspiral 
events will give a chirp mass spectrum 
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3. Determining the cosmic equation of state 



In order to explicitly show the dependence of the chirp 
mass spectrum on the cosmic equation of state, one 
needs to consider it in the context of general Friedmann- 
Robertson- Walker (FRW) cosmologies in the presence of 
a dark energy component. For general FRW cosmologies, 



Z.-H. Zhu et al.: Determining the cosmic equation of state using future gravitational wave detectors 



10 - 



(0.3,0) 

(0.3,-1/3) 

(0.3,-2/3) 



O 



cu 



5 - 








1.2 1.25 1.3 1.35 
Jt/M Q 



1.4 



(°m- "J 

(0.3,0) 

(0.3,-1/3) 

(0.3,-2/3) 




1.2 1.4 1.6 1.8 

Jt/M Q 



Fig. 1. The chirp mass spectrum of NS/NS binary inspirals for a fiat universe, (fl m , &x) — (0.3,0.7), with ui x = (SCDM 
model), —1/3, —2/3 and —1 (ACDM model). The detector depth is assumed to be ro — 200 Mpc (left panel) and ro = 455 Mpc 
(right panel) respectively. Solid lines are for the nonevolution case while dash-dot lines are for the evolutionary case from 
Schneider et al. (2000). 



the metric of space time is described by (in the c = 1 
unit) : 



ds 2 



-dt 2 + R 2 (t)[d X 2 + f 2 ( X )(de 2 
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where f(x) = X for a nat universe (k = 0), f(x) = sinx 
for a closed universe (k = +1), and f(x)= sinh% for an 
open universe (k = — 1). The comoving distance x is (Zhu 
1998, Zhu & Cao 1999) 
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for LIGOIlEJ (right panel), respectively. In our calcula- 
tions, we use a flat universe: (£l m ,£l x ) — (0.3,0.7) with a 
constant u> x e (—1,0), the dimensionless Hubble constant 
h = 0.65 and the signal-to-noise ratio threshold po = 8. 
As it is shown, there exists a sensitive dependence of the 
chirp mass spectrum of NS/NS binary inspiral on the cos- 
mic equation-of-state. In general the smaller cu Xl the more 
compressed the spectrum, and the smaller the tail at large 
M . This fact may provide an independent and robust con- 
straint on the cosmic equation of state. 
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where i includes each component of matter and energy in 
the universe, u>i = Pi/pi parameterizes the effective equa- 
tion of state for the zth-component, f2j = (87rG/3iJg)pi07 
and flk = —k/R^Hft with Hq being the Hubble constant 
(Hq = lOO/ikmsec^Mpc" 1 , 0.5 < h < 1) and zero sub- 
scripts referring to the present epoch. Various distance 
measures, including Dl and cdt/dz that are essential to 
calculate the chirp mass spectrum Eq.(||) can be easily in- 
ferred from the above comoving distance (Zhu 1998, Zhu 
& Cao 1999). 

Fig. |l| explicitly demonstrates how the chirp mass spec- 
trum depends sensitively on the equation of state lo x with 
both cases of NS/NS merger rate nonevolution (solid lines) 
and evolution (dash-dot lines). We have use two detector 
depths, r = 200 Mpc for LCGT (left panel) and 455 Mpc 



4. Discussion 

The solid lines of Fig. [j] are for the case of NS/NS merger 
rate nonevolution, F(z) = 1. Unfortunately the evolution 
rate is probably much strong. The dash-dot lines of Fig. [l] 
are based on the result of Schneider et al. (2000) (Eq.(|l])). 
As is shown, the merger rate evolution makes the spec- 
trum less compressed and extends the tail towards large 
M, which is exactly opposite to the effect for lowering 
u> x - Because uncertainties in the merger rate evolution di- 
rectly lead to uncertainties in determining the cosmologi- 
cal equation of state, the contraints on dark energy from 
gravitational wave detectors could not be robust until we 
understand well the cosmic evolution of the NS/NS binary 
merger rate. 



1 http://www.ligo.caltech.edu/ ligo2 /index. html 

http://www.ligo.caltech.edU/docs/T/T990080-00.pdi 



Z.-H. Zhu et al.: Determining the cosmic equation of state using future gravitational wave detectors 



10 - 



o 



cu 



5 - 







i — i — i — | — i — i — i — i — | — i — r~ 

(0.2,-1/3) 
(0.2,-2/3) 
(0.2,-1) 

(0.4,-1/3) 
(0.4,-2/3) 
(0.4,-1) 




Po= 8 
s r =200 Mpc 

^ O =1.2M 



1.2 1.25 1.3 1.35 1.4 
Jt/M Q 



(0.2,-1/3) . 

(0.2,-2/3) 

(0.2,-1) 

(0.4,-1/3) 
(0.4,-2/3) 
(0.4,-1) 




Po= 8 
r =455 Mpc 



Fig. 2. The chirp mass spectrum of NS/NS binary inspirals for a fiat universe with different matter density fi m = 0.2 (solid 
lines) and Q m — 0.4 (dash-dot lines). The detector depth is assumed to be ro = 200 Mpc (left panel) and ro = 455 Mpc (right 
panel) respectively. 



Another uncertainty of the method proposed here 
comes from the degenerate between the cosmic equation 
of state, u) x , and the cosmic matter density, fL„. Fig. g ex- 
plicitly shows how Q m changes the spectrum. A spectrum 
for a flat universe with lower n m and larger lo x may over- 
lap with the one for a flat universe with higher il m and 
smaller u x . One need an independent cosmological test to 
break down this degenerate, such as clusters of galaxies, 
gravitational lensing or the anisotropy of CMBR etc. 

Last but not least, neutron star masses are not all iden- 
tical (1.37M Q ), but distribute between a lower bound (e.g. 
1.29 M Q ) and a upper bound (e.g. 1.45 M Q ). This would 
introduce another probability function into Eqs.(3) & (4). 
Using a uniform distribution model, Finn (1996) shows 
that, the chirp mass spectrum will broaden symmetrically 
as the mass distribution broadens. Although the effect is 
different from the cosmic equation of state, the uncertainty 
of neutron star mass distribution will also lead the un- 
certainty into the constraints on dark energy from future 
gravitational wave events catalogue. 

As is shown above, the constraints on the cosmic equa- 
tion of state from the chirp mass spectrum of NS/NS 
merger events is promising, but it needs to appeal to 
statisitics, evolution and neutron star mass distributions. 
It is worth mentioning that, if NS/NS merger events are 
observed to be associated with gamma-ray bursts, one 
could measure the redshift of these events through the 
galaxies hosting the associated GRBs. Then a few NS/NS 
merger events would determine the cosmic equation of 
state as what the high-z SNIa has done. 
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